Signature of directed chaos in the conductance of a nanowire 
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We study the conductance of chaotic or disordered wires in a situation where equilibrium transport 
decomposes into biased diffusion and a counter-moving regular current. A possible realization is a 
semiconductor nanostructure with transversal magnetic field and suitably patterned surfaces. We 
find a non-trivial dependence of the conductance on the wire length, ft differs qualitatively from 
Ohm's law by the existence of a characteristic length scale and a finite saturation value. 
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According to Ohm's law, the resistance of a wire is 
proportional to its length. This is a straightforward con- 
sequence of the diffusive motion of electrons in the disor- 
dered potential of a normal material. However, unlike the 
time when Ohm arrived at his fundamental observation, 
conductors can be tailor-made today with almost com- 
plete control over the microscopic structure. Therefore 
it is important to understand the consequences of non- 
diffusive electron dynamics on the electronic conductance 
or other transport properties. This question has been 
studied in much detail for semiconductor nanostructures 
in which the motion of electrons is ballistic rather than 
diffusive 0,0>0>0- In such systems disorder is negligible 
and consequently all transport properties are determined 
by the shape of the sample, as in a billiard model. For ex- 
ample, in the ideal case of a perfectly clean nanowire with 
parallel walls the resistance should be zero independent 
of the length, and indeed this remarkable prediction has 
been confirmed experimentally Q. Beside ballistic sys- 
tems, also the effects of anomalous diffusion on the elec- 
tronic or thermal_conduction properties have attracted a 
lot of attention 



In the present paper we study the electronic conduc- 
tance of a wire in the case of a different and very profound 
modification of the microscopic dynamics. We consider 
systems where directed chaos leads to biased diffusion in 
the absence of any potential gradient. Directed chaos 
means that the time-averaged velocity of chaotic trajec- 
tories is non-zero due to broken time-reversal symmetry 
and due to the specific phase space structure. This effect 
may occur in various types of systems including, e.g., cold 
atoms in suitably pulsed optical potentials or chains of 
electronic billiards in a transversal magnetic field. It has 
been investigated both, theoreticallyand ex peri mentally, 
in a number of recent publications H, EJ EiEl H] ■ The 
interest is due to some intriguing and potentially very 
useful properties. For example, quantities like velocity 
average, velocity dispersion or scattering delay times are 
intrinsically dependent on the transport direction. How- 
ever, all previous studies focused on the ratchet-like di- 
rected transport in effectively infinite periodic systems 
with directed chaos. In contrast, we address here for the 



first time the typical electronic setup of a finite sample 
which is coupled to two electron reservoirs. We show 
that directed chaos has profound consequences also in 
this context where the transport velocity is not directly 
measurable; instead the conductance becomes the most 
basic and most relevant quantity. Due to biased diffu- 
sion, a new length scale A c h appears and rules the asymp- 
totic decay of the conductance with the sample length, 
see Eq. (|10|) below. Moreover, chaotic trajectories can 
propagate through samples of arbitrary length and conse- 
quently the conductance approaches a non-zero constant 
given in Eq. © . Note that it is not possible to reduce the 
description of directed chaos to a diffusion equation with 
bias. In our explicit result for the conductance, Eq. {T7J|, 
we must account also for the detailed structure of the 
underlying mixed phase space. 

To be specific, we investigate the prototypical model 
first introduced in 0. We consider a two-dimensional 
electron gas (2DEG) confined to a quasi one-dimensional 
channel (Fig. [TJi). One wall of the channel is straight. 
Electrons are specularly reflected, but no back scatter- 
ing occurs in the transport direction. The other wall has 
a rough surface causing strong and essentially random 
scattering. The details of this roughness are not crucial 
(see below). Directed chaos is induced by a perpendicu- 
lar magnetic field which breaks time-reversal invariance. 
We stress that a realization of our model does not require 
more than a novel combination of elements which are 
all well understood and experimentally accessible in the 
context of a mesoscopic 2DEG, namely transversal mag- 
netic fields of moderate strength, negligible bulk scatter- 
ing and surfaces which are either disordered or manufac- 
tured with a precisely defined geometry 0, 0, E| 0, H, E| . 

For our purpose it is sufficient to treat the electrons 
as independent classical particles, i.e., we assume a semi- 
classical regime with many transversal modes in the chan- 
nel, N ~ m*v-pb/h 3> 1. We use dimensionless units 
in which channel width and Fermi velocity are unity, 
b = v-p = 1. The strength of the magnetic field is para- 
meterized by the cyclotron radius r c = m*VY / eB. Besides 
the length of the sample this is the only free parameter in 
our problem. We will assume that the magnetic field is 
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FIG. 1: (a) A prototypical example for directed chaos is a 
2DEG in a perpendicular magnetic field which is confined by 
two parallel walls with different surfaces. The upper wall re- 
flects specularly while at the lower wall the reflection angle is 
essentially random. Possible physical realizations of this ran- 
domness are magnified in the insets (see text). Typical tra- 
jectories are either regular (dashed line) or random (full line) 
and transport in opposite directions. In (b) the transversal 
Poincare sections for v x > and v x < are displayed. 

not too strong, r c > 1. In this case there are no pinned 
orbits in the bulk of the channel and the phase space 
contains only two types of electron trajectories: regular 
orbits skipping along the clean channel boundary, and 
chaotic or random orbits which are reflected from both 
walls (Fig. ^). The regular orbits are transporting con- 
tinuously in one direction, say to the left. On the average, 
the chaotic orbits are transporting in the opposite direc- 
tion, thus compensating for the regular transport and 
making the system unbiased as a whole. However, the 
transport velocity x of a chaotic electron is fluctuating, 
in fact the dynamics is diffusive with a superimposed drift 
along the channel. The average drift velocity can be ob- 
tained by application of a phase space sum rule [||. For 
our model we find for the long-time velocity average of 
almost all chaotic trajectories 

1 9 — sin 8 cos 8 

Vch ~ 2 tt(1 - cos 0) - sin 8 + 8 cos 8 ^ ' 

with 8(r c ) = arccos(l — l/r c ) This result is inde- 

pendent of the precise modelling of the rough channel 
surface as long as the phase space structure of Fig. ^ is 
preserved. In 12] we considered the two extreme cases 
shown in the insets of Fig. QJl. In one case the surface 
is a periodic array of semicircular scatterers with small 
radius R — > i.e., the dynamics of the system is 
deterministic and there is no disorder whatsoever. In 
the second case the direction ip of the trajectory was 
randomized upon every scattering from the rough sur- 
face. Specifically it was chosen with probability density 
P(<p) — ^sin^j from the interval [0, 7t] such that the in- 
variant measure on the energy shell dxdydip was pre- 
served. A physical realization of this behavior is a dis- 
ordered surface with a correlation length that is below 



the Fermi wavelength. While this second system is non 
deterministic when approximated classically, its trans- 
port properties are essentially the same as for the case 
of deterministic directed chaos. In particular, in both 
cases an ensemble of chaotic trajectories spreads diffu- 
sively around the moving center of mass, (Ax 2 ) = D c ht 
[l2| . The precise value of the diffusion constant depends 
on the detailed modelling of the rough boundary. Ana- 
lytical results for D c h are available in the case of random 
scattering ^| and therefore we shall use this version of 
the model in the numerical calculations below. 

The electronic conductance is obtained within the 
framework of the Landauer-Buttiker formula G(L) — 
(2e 2 /h)T(L). We approximate the transmission semi- 
classically, T ~ Nt(L), and discuss quantum corrections 
at the end of this letter. t(L) denotes the total classi- 
cal probability that an electron is transmitted through a 
sample of length L if it enters the system at x = from 
the left with random initial conditions P(y,ip) = ^cos ip 
(y G [0,1], ip G [— §,+§])■ According to our convention 
(regular orbits are skipping to the left) these initial con- 
ditions are within the chaotic component of phase space. 
Similarly, we define the probability t'{L) that an electron 
is transmitted if it enters at x = L from the right with 
an analogous distribution (but tp e [? , ^-]). The total 
transmission probability must be the same for the two 
distinct transport directions, 

t(L) = t'(L) . (2) 

There are various ways to arrive at this fundamental iden- 
tity. For example, one observes that t ^ t' would result 
in the accumulation of particles in one of the reservoirs 
even when the system is in thermal equilibrium. Alterna- 
tively, a microscopic derivation can be based on the fact 
that the scattering map of Hamiltonian systems is area 
preserving. Although Eq. (J2J) is not specific for directed 
chaos, this identity has very interesting consequences in 
the present context. While t(L) is entirely due to chaotic 
trajectories whose properties are not immediately acces- 
sible, t'(L) can be decomposed into conditional probabil- 
ities for regular and chaotic trajectories, 

t'{L) = fi rcg t' Icg (L) + n ch t' ch (L) , (3) 

where 

If,, ie-sia8 cosd 

^ = 2 J Icg dydip C ° SIP = 2 l-cos0 (4) 

is the the relative area of the regular component in the 
transversal Poincare section (Fig.^>), and /x c h = 1 — Mrcg- 
Due to the lack of back scattering along the regular tra- 
jectories we have t[ eE (L) = 1. Moreover, it is clear that 
t' ch — > for L — > co. This is so because for trajecto- 
ries which are moving through arbitrarily long samples 
against the average chaotic flow the time-averaged ve- 
locity cannot converge to i> c h- Hence these trajectories 
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must be of measure zero in phase space. Taken together, 
the mentioned facts yield a remarkable result which is 
illustrated in the right inset of Fig. |2 for long systems 
the probability that a chaotic trajectory transmits from 
x = to x = L is given by the relative phase space area 
occupied by the counter-moving regular trajectories, 



t(L) — > jltreg (L — > oo) , 



(5) 



The goal is now to understand quantitatively how t' ch (L) 
decays from t' ch (Q) = 1 to zero as the length of the sys- 
tem increases. The results of numerical simulations for a 
wide range of values r c are shown in Fig. El The data sug- 
gest an exponential behavior for long wires. This can be 
understood after replacing the microscopic dynamics by 
the Fokker-Planck equation (FPE) for biased diffusion. 
At x = and x = L we assume absorbing boundary 
conditions. The probabilities to reach these exits from a 
point < x < L within the sample are then 



Pl{x) = 
respectively 0. Here 



-z/A c h 



-i/Ach 



=-L/X ah 



I _ g-^/Ach 
1 _ p-i/Ad, ' 



Ach = Dch/v, 



ch 



(6) 



(7) 



(8) 



denotes the Peclet length of the biased diffusion process. 

In order to contribute to the transmission from the 
right to the left end of a wire of length L, a chaotic par- 
ticle should first be transmitted through a segment of 
length I < L and then, starting from x — L — I, be ab- 
sorbed at x = 0. Based on this argument we propose as 
a recursion relation for the chaotic transmission 



t' ch (L)=t' ch (l)p' L (L-l). 



(9) 



Note that this relation cannot be valid for arbitrary I. 
For example, I = leads to t' ch (L) — as for a diffus- 
ing particle starting at one of the absorbing boundaries 
the probability to reach the opposite end is identically 
zero. The reason for this limitation can be understood 
as follows. Upon replacing the original dynamical sys- 
tem by a ID FPE we have discarded all information 
about the momentum of the electron. Strictly speak- 
ing it is then impossible even to define a transmission 
probability since this requires to enter the system with 
a given direction. Hence, for Eq. © to be valid, the 
length scale for momentum correlations should be neg- 
ligible compared to the distance from the boundaries, 
A c h x and A c h <C I. Under this assumption we find to 
leading order t' ch (L) = t' ch (l) e - ( i- ')/ Ach with the solution 



t' ch (L) = c exp(-L/A 



ch) 



(10) 



The dashed line in Fig. |2 shows this exponential for 
r c = 50 (with suitably chosen prefactor). Indeed the 



asymptotic decay of the transmission probability is re- 
produced. However, for short systems the behavior is 
clearly not exponential. Moreover, even in the asymp- 
totic regime L — ► oo the FPE approach is not accurate 
enough to predict the prefactor c of the exponential decay 
This is no surprise. According to Eq. i fTUj l. a differ- 
ent prefactor corresponds to an additive constant in the 
system length L and this type of error must be expected 
after discarding the correlation length of the momentum. 
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FIG. 2: The chaotic transmission (+) is compared to Eq. I16H 
for r c = 200, 50, 15, 10, 6, 3, 2 (top-bottom). Each data point 
represents 10 6 trajectories. For r c = 50 the asymptotic expo- 
nential is shown with a dashed line. Left inset: Dependence 
of the Peclet length on the cyclotron radius. Right inset: For 
r c = 10 the total transmission (•) is compared to Eq. 1171 



and to the asymptotic constant fj, r , 



0.294. 



We conclude that a satisfactory theory for t' ch (L) can- 
not be based on the FPE alone. As an alternative model 
let us consider a generalization of the persistent ran- 
dom walk 14] to a biased walk (BPRW): A particle 
moves with velocity ±wp and is reflected with direction- 
dependent probabilities at obstacles with spacing L\. For 
a segment of length L = nL\ we denote the transmission 
and reflection probabilities by t n = 1 — r n = t(nL\) (left 
to right) and t' n = 1 — r' n = t' ch (nLi) (right to left). 
A multiple-scattering expansion allows to express these 
probabilities in terms of a single segment. We find 

1-r'Jn 



t n — 



t = 



{t'JhY-r'Jn 
1-ri/rJ 



(11) 



(12) 



(*i/*i) n -n/ri ' 

In our case the chaotic transport is biased to the right 
Therefore we assume t\ > t[ and find in the limit n —> oc 



too = l - n/r 



and 



Ki — too (h/t'i)' 



(n 



oo 



(13) 



(14) 
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However, there is no direct connection between the pa- 
rameters of the BPRW and the dynamics of our orig- 
inal model. In order to close this gap we must make 
use of the information about the underlying phase space 
structure, Eq. and the result obtained within the 
FPE approach, Eq. (JTJJJ. The former implies t^ = fi lcg 
or equivalently ri/r[ = /i c h- Further the comparison of 
Eqs. (HU and (JTU) yields 



c exp(-n£i/A ch ) = Mrog 



Now it is easy to read off the correct prefactor of the 
asymptotic exponential decay, c = /irog, which was used 
for the dashed line in Fig. [21 On the other hand we infer 
exp(Li/A c h) = ti/t'i which is substituted into Eqs. i(TT|) 
and i)12fl . The final result for the chaotic transmission 
probability from right to left is then 



= 



Mreg 



exp(L/A ch ) - Mch 



while the total transmission is given by 
t(L) = 



1 - /i ch exp(-L/A ch ) 



(16) 



(17) 



This non-trivial prediction, which does not contain any 
free parameters, is confirmed numerically in Fig. The 
fact that an appropriate synthesis between FPE and 
BPRW should be used to reproduce the data was not 
at all obvious. Note that in an analogous approach to 
the unbiased case the PRW yields Ohm's law in the form 
t(L) = (D/v^L- 1 (L -> oo) while the FPE implies only 
t ~ L~ l . In contrast, for the biased case both approaches 
contribute complementary information and are in fact 
mutually incompatible approximations |l8| . 

Small but systematic deviations are visible in Fig. [2] 
for both, very small and very large cyclotron radius. At 
least partially these deviations can be attributed to di- 
rect trajectories which escape from the system before be- 
ing scattered. Such trajectories are always sensitive to 
the details of a given model and therefore not of primary 
interest here. We restrict the discussion to some sim- 
ple examples. Trajectories which enter at x = with 
a steep angle — ir/2 < ip will reach x = again after 
completing a simple arc. From the geometric condition 
cos ip < y/2r c we see that the corresponding phase space 
volume vanishes as r~ x . Therefore these trajectories lead 
to deviations for small cyclotron radius. For short sys- 
tems or for large cyclotron radius and shallow incidence 
ip ss there can also be directly transmitted trajectories. 
A rough estimate requires L < \/8r c for their existence 
which is compatible with Fig. [21 (L — 40 for r c — 200). 

According to Eq. I|17|) the conductance of a wire with 
directed chaos saturates to a finite value as L — > oo. 
This claim is in sharp contrast to quantum localization 
which leads to vanishing conductance in any coherent 
quasi ID quantum system with uncorrelated disorder. 



These two contradicting statements can be reconciled as 
follows. In the presence of directed chaos the localiza- 
tion length £ diverges exponentially with the number of 
transversal modes (ln£ ~ N ~ h -1 ) ^^|. Thus, even for 
moderate N, the localization length will easily exceed the 
sample length L or the coherence length of the given ma- 
terial. In this regime we can safely ignore localization. 
Other quantum effects like weak localization or tunnel- 
ing are not expected to change our results qualitatively 



(15) although they may lead to small corrections ~ N 



the relevant parameters A and /i rog . A numerical analysis 
of such effects will be attempted elsewhere. 
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